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ABSTRACT 

The power spectra of black hole binaries have been well studied for decades, giving a 
very detailed phenomenological picture of the variability properties and their correlation with 
the energy spectrum (spectral state) of the source. Here we take the truncated disc/hot inner 
flow picture which can describe the spectral changes, and show that propagating mass accre- 
tion rate fluctuations in the hot flow can match the broad band power spectral properties seen 
in black hole binaries, i.e. give approximately band limited noise between a low and high fre- 
quency break. The low frequency break marks the viscous timescale at the outer edge of the 
hot inner flow, which is the inner edge of the truncated disc. The fluctuations in mass accretion 
rate propagate towards the central object in a finite time meaning the high frequency break 
is more complex than simply the viscous timescale at the inner edge of the hot flow because 
fluctuations on timescales shorter than the propagation time are incoherent. The model also 
predicts the Lense-Thirring precession timescale of the hot flow, as this is set by the com- 
bination of inner and outer radius of the flow, together with its surface density which is self 
consistently calculated from the propagating fluctuations. We show that this naturally gives 
the observed relation between the low frequency break and QPO frequency as the outer radius 
of the flow moves inwards, and that this model predicts many of the observed QPO properties 
such as correlation of coherence with frequency, and of the recently discovered correlation of 
frequency with flux on short timescales. 

We fit this total model of the variability to a sequence of 5 observed power spectra from 
the bright black hole binary XTE J 1550-564 as the source transitioned from a low/hard to 
very high state. This is the first time that a power spectrum from a black hole binary has been 
fit with a physical model for the variability. The data are well fit if the inner radius of the flow 
remains constant, while the outer radius sweeps inwards from ^ 75 — l2Rg. This range of 
radii is the same range as required by models of the energy spectral evolution, giving the first 
self consistent description of the evolution of both the spectrum and variability of BHB. 
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1 INTRODUCTION 

Emission observed from Black Hole Binaries (BHBs) is variable on 
both long and short timescales. On the longest timescales, changes 
in mass accretion rate drive changes in spectral state. At the low- 
est luminosities, the system is typically seen in the low/hard state, 
where the spectrum is dominated by a hard (photon index F < 2) 
power law tail, peaking at 100 keV. As the source brightens, there 
is an increasing contribution from the disc at low energies and the 
tail softens towards F ~ 2 (intermediate state). The tail can then 
either remain strong as the source flux increases, giving the very 
high state (F ~ 2.5) or it can carry only a very small fraction of the 
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luminosity, giving the high/soft state (F ^ 2.2) (e.g. see reviews by 
Remillard & McClintock 2006; Done, Gierlinski & Kubota 2007, 
hereafter DGK07). 

Much of this long term spectral evolution can be successfully 
described by the truncated disc model (DGK07). Here, the assump- 
tion is that the cool, geometrically thin, optically thick disc is trun- 
cated at some radius, ro, which is much larger than the last sta- 
ble orbit, riso at the lowest luminosities. The inner accretion flow 
extends from Vo to an innermost radius, r^, forming some sort of 
hot, geometrically thick, optically thin (optical depth r ~ 1.5) 
flow, similar to an Advection Dominated Accretion Flow (ADAF, 
Narayan & Yi 1995). Decreasing the truncation radius of the thin 
disc with increasing luminosity (and therefore average mass accre- 
tion rate) leads to a stronger disc component and greater illumina- 
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tion of the hot flow by the cool disc photons and hence to a softer 
Comptonised tail. This spectral evolution has a natural end when 
To ~ Tiso, marking the transition to the disc dominated soft states 
seen at higher luminosities. 

There is also rapid variability on timescales of 100-0. Is 
which is almost exclusively linked to the Comptonised tail (e.g. 
Churazov et al 2001). This is as expected as the disc emis- 
sion can only vary on a viscous timescale of t^isc = 4.5 x 
lO~^a'^{h/r)-^{R/&Rgf'^{M/lQMQ) s i.e. 500 s at the last 
stable orbit for fiducial thin disc parameters of h/r ~ 0.01 and 
a = 0.1. 

A power spectral analysis of the rapid variability shows that 
it is generally composed of a broadband noise continuum, with a 
low frequency quasi-periodic oscillation (LF QPO) superimposed 
on this. The broadband noise can be very roughly characterized as a 
twice broken power law, with fPf oc / at the lowest frequencies, 
breaking to fPf oc /" above the low frequency break /;,, and then 
breaking again to fPf oc /^^ above a high frequency break at /;,. 
The low frequency break is tightly correlated with the spectral state 
of the source, with ft ~ 0.01 Hz for the dimmest low/hard states, 
and increasing to ~ 1 Hz as the source softens through the bright 
low/hard states into the intermediate state. The LF QPO properties 
are also tightly correlated with this power spectral and energy spec- 
tral evolution, with QPO frequency fgpo ~ 10/b increasing from 
0.1-10 Hz, as the coherence, power and harmonic content of this 
feature also increase (e.g. van der Klis et al 2004). Again, this can 
be interpreted qualitatively in the truncated disc picture if the low 
frequency break and QPO are somehow associated with the trunca- 
tion radius. As the disc extends further inwards, the spectrum soft- 
ens, and all frequencies associated with this radius increase (e.g. 
DGK07). 

However, to go beyond this into a quantitative description re- 
quires a specific model for both the broadband noise and QPO. This 
is problematic, despite both having been known for many decades 
(e.g. van der Klis 1989). There are multiple potential models for the 
LF QPO in the literature which fall into 2 main categories: those 
associated with a geometrical misalignment of the accretion flow 
and black hole spin (Stella & Vietri 1998; Fragile, Mathews & Wil- 
son 2001; Schnittman 2005; Schnittman et al 2006; Ingram, Done 
& Fragile 2009, hereafter IDF09), and those associated with wave 
modes of the accretion flow (Wagoner et al 2001; Titarchuk & Os- 
cherovich 1999; Cabanac et al 2008). Most of these concentrate on 
matching the QPO frequency, but the spectrum of the QPO gives 
additional constraints. This is similar to that of the spectrum of the 
broadband variability, showing that they both arise predominantly 
from the Comptonising region rather than the disc (e.g Gilfanov 
et al 2003; Sobolewska & Zycki 2006), favoring models in which 
the modulation arises directly from the Comptonised emission e.g. 
IDF09, where the QPO is set by Lense-Thirring (vertical) preces- 
sion of the entire hot inner flow interior to the disc truncation radius 
at r-Q, or by a mode of the hot inner flow (Cabernac et al 2008). 

The physical origin for the viscosity of the flow is the 
Magneto-Rotational Instability (MRI: Balbus & Hawley 1991). 
This is inherently variable, with large fluctuations in all quantities, 
both spatially and temporally (Krolik & Hawley 2002), making it a 
natural origin for the broad band noise (Noble & Krolik 2009; In- 
gram & Done 2010; Dexter & Fragile 2011). However, these fluc- 
tuations also effectively shred any coherent wave modes in the flow 
(see e.g. Reynolds & Miller 2009; Henisey et al 2009). This effec- 
tively rules out trapped wave propagation as the origin of the LF 
QPO, leaving Lense-Thirring precession as the most likely candi- 
date. Precession of the entire hot flow from ro to Vi can match the 



observed LF QPO frequency in both BHBs (IDF09) and neutron 
stars (NS; Ingram & Done 2010), and provides a clear mechanism 
to match the spectrum as this is a modulation of the Comptonising 
region. 

Thus the entire power spectrum can be explained by MRI fluc- 
tuations in a hot flow, which is also precessing around the black 
hole. However, the power spectrum does not represent all the in- 
formation contained in the variability, as it uses only Fourier am- 
plitudes, not phases. This is important as the light curves contain 
additional correlations which give a linear rms-flux relation (Utt- 
ley & McHardy 2001). Short segments of a longer light curve have 
a mean, /, and variance, a^, which are related such that I oc a 
(after binning: Uttley & McHardy 2001). This is equivalent to the 
flux on these timescales having a log-normal distribution (Negoro 
et al 2000), and rules out simple models of the variability where 
the light curve is made from adding together multiple, uncorrelated 
events (Uttley & Mchardy 2001; see also DGK07). Instead, this 
can be produced if the light curve is made from a multiplicative 
process, rather than an additive one. Again, the MRI in the hot flow 
gives a physical interpretation to this. The MRI at large radii pro- 
duces intrinsic fluctuations in the density of the flow. These fluctua- 
tions propagate down to smaller radii on a viscous timescale, so all 
higher frequency fluctuations are smoothed out. These smoothed 
and lagged fluctuations in mass accretion rate modulate the MRI 
fluctuations produced by the next radius, and so on, down to the 
smallest radii in the flow (Luybarski 1997: hereafter L97). This 
naturally produces a light curve which has an rms-flux relation/log- 
normal flux distribution (Kotov et al 2001: hereafter KOI; Arevalo 
& Uttley 2006: hereafter AU06; also see Misra & Zdziarski 2008). 

Here we use these ideas to build a model for the entire power 
spectrum, where the broadband noise arises from propagation of 
MRI fluctuations through the hot flow from ro to and the LF 
QPO arises from Lense-Thirring precession of the same hot flow. 
We fit this model to a sequence of power spectra of XTE J1550-564 
observed during a transition from low/hard to intermediate state, 
and show that these give a good quantitative fit for Vo moving from 
~ 75 — 10i?g, as also required to explain the sequence in energy 
spectra during this transition. This gives very strong support to the 
underlying truncated disc geometry, and for Lense-Thirring preces- 
sion as the physical origin for the LF QPO. 



2 THE SIMPLIFIED MODEL 

We first introduce the simplified toy model considered in IDF09 
and Ingram & Done (2010), whereby the low frequency break, 
fb, occurs at the viscous frequency of the truncation radius and 
the QPO frequency, fgpo, is the precession frequency of the 
flow. Using the Shakura Sunyaev (1973) viscosity prescription, 
fb = fvisc{ro) = aih/r)"^ fk{ro) = -Vr{ro)/Ro where a is 
the viscosity parameter, h/r the flow semi-thickness, fk the Kep- 
lerian frequency and v,. is the infall velocity. The Lense-Thirring 
precession frequency is then a weighted average of the point parti- 
cle precession frequency at each radius in the hot flow, so 

/;; fLTfk^r'dr 

fprec — fgPO — V'-o t ^ 'iJ 

J^^ fk^r^dr 

(Liu & Melia 2002) where Vi is the innermost point of the flow (i.e. 
the surface density is negligible interior to this), E is the surface 
density and 
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Figure 1. The QPO-break relation plotted in dimensions of c/Rg for a 
fiducial mass of 10 and IAMq for black holes (black circles) and neutron 
stars (red triangles) respectively. The fact that these frequencies lie on the 
same relationship for the two objects implies a common physical origin. 
The lines are predictions of the simplified model for spins of a, = 0.2 
(dot-dashed), 0.5 (dashed) and 0.998 (solid) with To ranging from to 
100 and = rj^,. For the blue lines, we assume the viscous frequency 
to be proportional to the Keplerian frequency. For the green lines, we as- 
sume f^isc = Br~™' ff^ where B = 0.03 and m = 0.5 and find good 
agreement with the trend in the data. 



predicts frequencies which are fairly close to the observations, it is 
clear that the gradient of this model in log space is different from 
that observed. 

Plainly the assumptions above are very simplistic. Global an- 
alytic models of the hot flow with a standard ct viscosity do not 
have fvisc oc /k as they depart from the self-similar solutions at 
r < 100 due to the requirement that the flow becomes supersonic 
(Narayan, Kato & Honma 1997; Gammie & Popham 1998). Full 
numerical simulations also show that a is not constant (e.g. Fragile 
et al 2007; 2009). Ingram & Done (2010) also suggest that ( can 
change in neutron stars as the material piles up onto a boundary 
layer. However, the similarity between the mass scaled frequencies 
seen in neutron stars and black holes shown in Fig. [T] make this 
now seem unlikely to be an important effect as it would not affect 
the black holes. 

Here then we simply assume that a(/i/r)^ is a power law 
function of radius, so that fvisc = Br~"^fK- We choose values 
for B and m which allow us to match the data in Figure [T] We see 
good agreement with the observations for B = 0.03 and m — 0.5 
(green lines), again for a, — 0.2 (dot-dashed), a, — 0.5 (dashed) 
and a, = 0.998 (solid). We use this specific prescription for the 
viscous frequency in the following section. 



fhT = fk 



4a, 

„3/2 



+ 



3a« 



(2) 



is the point particle Lense-Thirring precession frequency for a di- 
mensionless spin parameter a* (Merloni et al 1999). Here, r is di- 
mensionless, expressed in units of Rg — GM/(?. Solving this 
assuming a power law form for surface density in the hot flow 
E cx r"'' between an inner and outer radius for the hot flow Vi 
and To (Fragile et al 2007, IDF09) gives 



Iqpo 



(5-- 20 a,[l-{ri/roY'^+^] 



^(1 + 20 5/2-C l/2+Cri 



(r,/r„)5/2-C] Rg 



(3) 



Hence this model predicts the relation between /q po and /t 
which can be compared to the multiple observations of these fre- 
quencies in both black holes and neutron stars (e.g. Wijnands & van 
der Klis 1999; Klein- Wolt & van der Klis 2008). The observed rela- 
tion is continuous, implying that these frequencies show the same 
behavior in both sources, i.e. that neither can depend strongly on 
any property of the neutron star surface but are instead set by the 
accretion flow itself. We re-plot this data in Figure [T] normalizing 
the frequencies by mass for a fiducial mass of 1.4 and IOMq for 
neutron stars (red triangles) and black holes (black circles), respec- 
tively. This shows even more clearly that the two different types of 
object show the same observed relation between these frequencies 
as they now occupy the same range. 

The blue lines show the prediction of the toy model, where 
the hot flow has constant a — 0.2 and h/r — 0.2, surface density 
constant with radius (i.e. C = 0: Fragile et al 2007) between Vo and 
Ti, where is given by the bending wave radius. Warps in a large 
scale height flow are communicated via bending waves which have 
wavelength A cx r^^^ and so are smooth at large r and oscillatory at 
small r. The bending wave radius (r, = 3.0(/i/f')"^''^a*''^ Frag- 
ile et al 2007; 2009; Fragile 2009; IDF09) marks the transition be- 
tween the two regimes. We show a* — 0.2 (dot-dashed), a* — 0.5 
(dashed) and a, = 0.998 (solid). While this very simple model 



3 THE FULL MODEL 

We consider a model where local fluctuations in the mass accretion 
rate of the flow propagate down towards the central object (e.g. 
L97; KOl). Our method mainly follows that of AU06, with a few 
small differences. 

We split the flow into annuli, characterized by a radius r-„ and 
width drn, with logarithmic spacing so drn/rn is a constant for all 



annuli from To to r^. As fvisc oc r 



-(m+3/2) 



in our prescription. 



constant drn/rn also implies constant df / f which AU06 show is 
the feature required to produce a linear sigma-flux relation. 

MRI fluctuations throughout the flow generate variability in 
all quantities. Fluctuations in mass accretion rate will most likely 
be damped by the response of the flow on timescales shorter than 
the local viscous timescale. We therefore assume that the gener- 
ated power spectrum of mass accretion rate fluctuations at radius 
r„ is given by a zero centred Lorentzian cutting off at the viscous 
frequency 



1 + (///.>.c(r„))2 



(4) 



where fvisc = ^^'"r/r — Br~"^fK as discussed at the end of 
section[2]and a tilde denotes a Fourier transform. 

We start at the outermost annulus, so ri — ro, and generate 
the time dependent fluctuations in mass accretion rate, Th{ri,t), 
from equation |4] using the method of Timmer & Koenig (1995). 
We normalize each m(r„,f) to have a mean of unity and frac- 
tional variability a/I = FvarV Ndec where Fvar is the fractional 
variability per decade in radial extent and Ndec is the number of 
annuli per decade in radial extent. Thus the mass accretion rate 
across the first annulus is M{r\,t) = Morn(ri,t) where Mo 
is the mean mass accretion rate. This then propagates inward to 
the second annulus, travelling a distance dr\, which takes a time 
tiag = Rg dr\/vr{r\). When it arrives at r2, it has been filtered 
by the response of the flow which we take from Psaltis & Norman 
(2000) to get 
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^l+[{dr„/r„)if/f,isc{r„W 



The mass accretion rate at the n* annulus is then given by 



(5) 



(6) 



M{rn,t) = Mf{rn-l,t ~tlag)m{rn,t) 

where tiag = Rg drn/vrirn). However, equation|5]only filters out 
fluctuations on much shorter timescales (by a factor dr/r) than the 
typical timescales generated in the annulus (equation |4j and so we 

can say Mf{rn,f) ~ M(rn,f ) to a very good approximation. 
The mass accretion rate at the n*** annulus is therefore given by 



M{r„,t) = M{r„^i,t-tiag)m{rn,t), 



(7) 



until the A*'' ^ annulus which is r^. 

To transform this into a light curve requires an emissiv- 
ity, e{r) such that the luminosity from each annulus is given by 
dL{r,t) = 1/2 M{r„,t)e{r)r„dr„c^ where we assume the emis- 
sivity e(r) oc r ''b{r) where b{r) is a boundary condition. For 
a Newtonian thin disc, 7 = 3 and we have the stress free inner 
boundary b{r) = 3(1 — \/ Tn/ri) but we note that the large scale 
magnetic fields present in the large scale height flow can give a 
stressed inner boundary condition b(r) — 1 (Agol & Krolik 2000; 
Beckwith Hawley & Krolik 2008). We also allow 7 to be a free 
parameter as the emission need not exactly follow the radial de- 
pendence of gravitational energy release as long as the total energy 
release is gravitational. A different emissivity for different energy 
bands gives a way for the model to predict frequency dependent 
time lags between hard and soft X-ray bands (KOI; AU06). This 
will be the subject of a subsequent paper. 



4 THE FIDUCIAL MODEL 

Figure [2^ shows 20 s of the resulting light curve for a fiducial set 
of input parameters for a black hole mass of M = IOMq and a 
spin of a* = 0.5. We assume — 2.5, Vo — 20, Fvar ~ 0.4, 
B = 0.03, m — 0.5, 7 = 4.5 with a stressed inner boundary 
condition (see Section [2}. We calculate the light curve with 2^^ 
time points, corresponding to ~ 4096 s (a typical length for an 
RXTE observation) of data on a time binning of 0.00097 s, and 30 
radial bins. 

Figure I2J3 (red) shows the PSD of this light curve, while the 
black and green points show the effect of changing the number of 
radial bins to 10 and 100, respectively. Clearly, the high frequency 
power is not well resolved with only 10 radial bins, while the dif- 
ference between 30 and 100 is very small. Hence in all that follows, 
we use 30 radial bins for each simulation. 

The PSD shows the same characteristic broadband noise fea- 
tures as are seen in the power spectra of black hole binaries, namely 
band limited noise, with low and high frequency breaks, peaking 
between 0.1-10 Hz. Figure[2l: shows the rms-flux relation for the 
fiducial light curve, derived from splitting this into 4s segments. 
As with the data (Uttley & McHardy 2001), we see a large scat- 
ter before binning (gray points) but, after binning (red crosses), we 
retrieve a linear flux-rms relation (AU06). 



5 THE TRUNCATED DISC/HOT INNER FLOW MODEL 

The major prediction of the truncated disc/hot inner flow model is 
that the spectral softening as the source brightens from a low/hard 
through to intermediate states is caused by the truncation radius 
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Figure 2. Top (a): A 20 second section of the simulated light curve. Mid- 
dle (h): The simulated power spectral density calculated using 10 (black), 
30 (red) and 100 (green) radial bins. We see that, 10 bins is not enough to 
resolve the high frequency power but 30 bins is a good approximation. Bot- 
tom (c): The rms-flux relation for the light curve shown. We see that this is 
linear as observed. 



of the thin disc moving inwards (e.g. DGK07; Gierlinski, Done & 
Page 2008). This radius also sets the outer edge of the hot flow, so 
this predicts that decreases also. 

Figure[3] shows the predicted PSD for Vo — 50, 20 and 10, as 
required to match the energy spectral evolution (and low frequency 
QPO: IDF09), with all other parameters held constant at the fiducial 
model values described above. The model predicts that decreasing 
the outer radius of the hot flow leads to less low frequency power, 
while the high frequency power remains constant. This is precisely 
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Figure 3. PSD calculated using the fiducial parameters with To = 50 
(black), 20 (red) and 10 (green) with total fractional variability generated 
per decade in radius, Fvar, held constant. This has the same characteristics 
as the observed PSD of the data as the source softens from a low/hard to in- 
termediate state, namely that the low frequency power drops while the high 
frequency power remains constant. 

what is seen in the PSD of the data (DGK07; Gierlinski, Nikolajuk 
& Czerny 2008). 

This is the first physical model of the power spectral behavior 
which naturally reproduces the observations. The low frequency 
break is close to the frequency of the viscous timescale at Vo, as 
proposed by e.g. by Churazov et al (2001); Gilfanov & Areief 
(2005); DGK07; Ingram & Done (2010). However, the high fre- 
quency break is not at the viscous frequency at ri. We explore the 
origin of the high frequency break below. 

5.1 Effect of propagation on the PSD shape 

Figure |4ji shows a model where the variability at each radius is 
a Lorentzian at the local viscous frequency (see equation |4j, but 
with no propagation so there is no causal connection between an- 
nuli. We show the PSD of the resulting m(r,i, t) functions from 5 
of the 30 individual radial annuli, from Vo to as the gray lines 
on Figure|4^. These peak, as expected, at fvisciro) and fvisc{ri). 
The total variability (red) is an emissivity weighted sum of these 
fluctuations, but since they are uncorrelated, the effect of this is to 
strongly dilute the total variability seen. This total PSD does have 
ff^ fvisc(ri) ~ 12Hz as our emissivity weighting strongly fa- 
vors the smallest radii, but > fvisciro) (~ lOHz and ~ 0.3Hz 
respectively). In fact, to achieve fPf oc /" as observed, we would 
have to assume a completely flat emissivity profile, which seems 
very unlikely. More fundamentally, such uncorrelated fluctuations 
cannot reproduce a linear rms-flux relation. 

This is in sharp contrast to a model where fluctuations propa- 
gate down in radius. The resulting PSD from the same set of radii 
are shown in Figure|4j3, where the power in each annulus increases 
strongly with radius as the MRI power generated in each annulus is 
modulated by the propagating fluctuations from all radii prior to it. 
The red line shows the resulting emissivity weighted power spec- 
trum from the total flow assuming that time lags between radii are 
negligible. This preserves the maximum correlation between vari- 
ability at different annuli i.e. gives the least dilution between fluctu- 
ations in different annuli. This is very different to that in Figure|4^, 
both in normalization and shape. The normalization is dramatically 
enhanced because the long timescale fluctuations are correlated to- 
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Figure 4. Top (a): The gray lines are the power spectra of 5 simulated 
rn{r„, t) functions. We simulate 30 of these functions but, for clarity, only 
plot 5 without showing the errors. The red line is the PSD of the light curve 
created by assuming there to be no propagation (i.e.M(r„, t) = rn{rn,t)) 
and an emissivity index of 4.5. Because the functions we sum over are un- 
correlated, the PSD of the light curve looks like the (weighted) sum of the 
100 PSDs with the only difference being the normalization. Bottom (b): The 
gray lines are now the power spectra of M{r„, t) functions, i.e. we now al- 
low propagation. These are correlated at low frequencies but not at high 
frequencies allowing the model to reproduce the observed linear sigma-flux 
relation. The red line is the PSD of the light curve if we do not consider the 
propagation time between annuli (i.e. t;„g = 0) and the blue line results if 
we do consider the propagation time. We see that the red line differs from 
the top plot in 2 ways: the normaUzation is much higher and high frequency 
noise is lost. However, much more high frequency noise is lost for the blue 
line indicating that considering lags reduces high frequency noise. These 
plots illustrate that the prediction from shot noise models such as the top 
plot that the observed high frequency break is the viscous frequency at the 
inner radius breaks down once we consider a more advanced model capable 
of reproducing other observational properties. 



gether, so at low frequencies the power from different radii add to- 
gether as they are in phase. This gives ~ fvisciro) ~ 0.3Hz as 
the correlated variability weighting to larger radii is stronger than 
the emissivity weighting to smaller radii. However, at the fastest 
timescales, the power is mainly generated at the smallest radii, so 
it does not correlate with any other fluctuations generated at larger 
radii, so is not enhanced in the same way. 

The blue line shows how time delays dramatically change the 
high frequency break as the propagation time prevents the mass 
accretion rate from two consecutive annuli from being correlated 
on time scales shorter than tiag- This reduces the correlation be- 
tween the fastest timescale variability, strongly suppressing high 
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Figure 5. The PSD calculated assuming b{r) = 1 and 7 = 4.5 (black), 
6(r) =stress free and 7 = 4.5 (red), b{r) =stress free and 7 = 3 (green), 
all with Ti = 2.5. The blue points are for b{r) =stress free, 7 = 3 and 
ri = 6. This illustrates that we can reduce the predicted high frequency 
noise by changing boundary condition, emissivity index or inner radius. 



frequency power. Thus in the propagating fluctuation model, the 
low frequency break is ft ^ fvisciro) but fh « fvisciri). 



5.2 Emissivity and boundary condition 

We use an emissivity to translate the fluctuations in mass accretion 
rate to a luminosity. This emissivity is in two parts, firstly a power 
law dependence in radius, and secondly a boundary condition. Our 
fiducial model parameters have 7 — 4.5 and a stressed boundary 
condition, fe(r) — 1. This emissivity peaks at ri, so fluctuations 
from the very smallest radii are given most weight. 

Figure[5]compares this (black line) with results using the same 
power law radial dependence, but with a stress-free inner boundary 
condition (red line), b{r) — 3(1 — \Jt jri). This emissivity goes to 
zero at the innermost radius, so the highest frequency fluctuations 
are strongly suppressed. However, this also has a more subtle effect 
on the region between the two breaks, as there is a gradual decrease 
in weighting of fluctuations below r = 2ri , and a stronger weight- 
ing to the fluctuations at larger radii, giving the tilt between ft and 

This effect is similar to that of changing the radial dependence 
of the emissivity. The green line shows 7 = 3 with a stress-free 
boundary condition, showing an even stronger tilt to the PSD be- 
tween ft and fh (green). However, it is also similar to changing 
the inner radius of the flow. The blue line in Figure[5]shows the re- 
sulting PSD from 7 = 3 and a stress-free inner boundary condition 
with ri = 6. Thus there are degeneracies between the two parts to 
the emissivity and the inner radius, making it unlikely that they can 
all be uniquely constrained by the observed PSD. 

It is clear from this analysis that while the low frequency break 
is fairly strongly linked to the viscous timescale of the outer ra- 
dius of the hot flow (as assumed in section [2l(, the high frequency 
break is rather more complex, depending on propagation correla- 
tions, emissivity, boundary condition and inner radius in addition 
to the viscous timescale. This makes it difficult to directly asso- 
ciate the high frequency break with any physical parameter of the 
models. Instead, we now use the additional information from the 
QPO to remove some of the degeneracies inherent in tiiis model for 
the broadband noise. 



6 THE QPO: PRECESSION AND SURFACE DENSITY 

For our fiducial model, we used the observed relation between the 
low frequency break and LF QPO to set the radial dependence of 
the viscous timescale, assuming that the low frequency break was 
set by the viscous timescale at Vo and that the QPO was Lense- 
Thirring precession of the entire hot flow from ro to ri (Section|2ll. 
This assumed that the surface density of the hot flow, E = Sor"'' 
between ro and ri, with ( = 0. However, the broadband noise 
model described above calculates a self-consistent surface density 
as mass conservation implies 



M{r„,t) = -27rr„«r(r„)E(r„,f), 



(8) 



(Frank, King & Raine 1992) where Vr is now expressed in units of 
c, M in units of Afo, E in units of Mo/(cRg) and r in units of Rg. 
Using our velocity prescription, we can then easily show 



E(r„,t) = 



B 



(9) 



This means that, for the time averaged surface density, ( — m— 1 /2 
giving extra physical motivation for the parameters used in section 
|2](C = Oandm = 1/2). 

Figure |6] shows E(rvi) plotted at a number of different times 
(0, 256, 512, 1792 s) along with the corresponding M(r,t) 
function. M{r,t) is quite clearly more variable at small r. This 
is because we have assumed the variability generated in each log- 
arithmic annulus to be the same but annuli at smaller radii include 
also the fluctuations that have propagated down from large r and so 
the emitted variability is greater (see Figure |4j5). We do not see a 
drop off in surface density at the bending wave radius like that seen 
in simulations (e.g. Fragile 2009) because we assume that the infall 
velocity can be given by a power law. It is clear that, for the surface 
density drop off at a given point, the infall velocity must accelerate 
at that point. In a future paper, we will investigate this model with 
a more advanced velocity prescription. 

Therefore the broadband noise model above, set by ri, ro and 
Fvar predicts the QPO frequency at any point in time. The fluctu- 
ations in surface density with time predict that the QPO frequency 
changes, i.e. it is quasi-periodic rather than truly periodic. How- 
ever, the precession frequency will not respond instantaneously to 
these changes, as their effect is only communicated across the en- 
tire hot flow by bending waves. These travel at the sound speed, 
faster by a factor ~ a than the viscous timescale across the region, 
so we calculate the QPO frequency every ~ 4 s rather than at every 
point. We then average these values to get the predicted QPO fre- 
quency, fqpo, and use the dispersion around this to set the r.m.s. 
variance of these QPO frequencies, ag po ■ 

Fig.|7]shows OQpo / Jqpo as a function of fgpo as ro varies 
from 300-10 in the fiducial model. This decrease in ro not only 
leads to an increase in QPO frequency, but also to a decrease in the 
QPO width, or equivalently, an increase in its coherence/quality 
factor Q — fqpo /o'qpo- This correlation is well known in QPO 
data from BHB (e.g. Belloni, Psaltis & van der Klis 2002; Rao et al 
2010). Our model provides the first physical explanation for this ef- 
fect as the smaller the radial extent, the higher the QPO frequency, 
but also the smaller the fluctuation power, giving smaller jitter in 
frequency. The red squares in Fig [7] show the ohsen'ed frequency 
and width of the QPO from data from the 1998 rise to outburst 
of the BHB XTE J1550-564 (see section 7). The model matches 
the trend in the data fairly well, and forms a lower limit to the 
width of the QPO. However, other effects such as the on-time of 
the QPO (see Lachowicz & Done 2010) can decrease the coher- 
ence of the signal, so here we simply use the data to determine the 
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Figure 6. Top: Mass accretion rate as a function of radius sliown here at a 
number of different times. Bottom: Surface density as a function of radius 
shown at the same times as the lines of corresponding colour in the top plot. 
This is calculated by applying mass conservation in the flow. The dotted 
line in the bottom plot represents the bending wave radius for a, =0.5 
and h/r = 0.2 which seems to trace the point where the surface density 
begins to drop off. 

QPO width (parametrized as a quality factor Q = fgpo /o'qpo) 
after using the broadband noise surface density to determine the 
QPO frequency, fgpo- 

The model also predicts another correlation, one between the 
QPO frequency and flux on short timescales. The top plot of Figure 
[8] shows this for the fiducial model (i.e. ro = 20), with precession 
frequency calculated every 4s together with the instantaneous lumi- 
nosity at that time. After binning (red crosses), there is a clear lin- 
ear relation between the two. This happens because both the QPO 
frequency and the luminosity depend on the mass accretion rate 
fluctuations. A perturbation in mass accretion rate at large r will 
lead to a perturbation in the surface density. This will reduce the 
precession frequency but will have little effect on the luminosity 
because the emissivity is quite steeply weighted towards small r. 
Later on, this perturbation will have propagated inwards to small 
r where it has the effect of increasing the precession frequency, 
but now also has much more of an effect on the luminosity. Heil 
et al (2011) have recently discovered this correlation in data from 
the 1998 rise to outburst of XTE J1550-564 (the same data as we 
will be considering in section [TJ. They also find that the gradient 
of this relation is steeper for observations with a higher QPO fre- 
quency. This is also predicted by the model as illustrated in the bot- 
tom plot of Figure [8] where we have measured the gradient of the 



Figure 7. Fractional variability of the precession frequency plotted against 
the average precession frequency (black line). These are calculated by mea- 
suring the precession frequency every 4s for a number of different trun- 
cation radii, ranging from 300 — 10, and taking the average and standard 
deviation over a 2048s duration. The red squares show the observed QPO 
width and frequency in data from the 1998 rise to outburst of XTE 1550- 
564. We see broad agreement with the data, however, other effects such as 
on-time of the QPO can decrease the coherence of the signal so we note that 
we are only able to predict a lower limit for the width of the QPO 



fgpo-L relation and the average QPO frequency for 11 different 
To values. There is clearly a very strong correlation as is seen in the 
data. This is because an absolute change in precession frequency 
depends on a. fractional change in mass accretion rate whereas an 
absolute change in luminosity depends on an absolute change in 
mass accretion rate. The same absolute change in mass accretion 
rate at a given radius and time constitutes a larger fractional change 
for small ro than for high ro - Therefore the luminosity will experi- 
ence exactly the same change in both instances but the precession 
frequency will undergo a larger change when Vo is smaller. The fact 
that these are predicted properties of the model constitutes strong 
support for its validity. 

In the following section, we will fit the PSD from our model 
to data. For this, we will need to predict a shape for the QPO light 
curve as well as a frequency. The data show that the QPO has a 
power spectrum which can be represented by a Lorentzian at the 
fundamental frequency, fqpo, together with its second and third 
harmonic and sub-harmonic i.e. at 2 fQpo, Sfgpo and 1/2/qpo 
(e.g. Belloni, Psaltis & van der Klis 2002). Our model for the QPO 
in terms of Lense-Thirring precession predicts the shape of the 
modulation of the emission from the hot flow via variation of pro- 
jected area, self-occultation and seed photons (IDF09). We will ex- 
plore this further in a later paper (Ingram, Done & Zycki in prep), 
but here we simply assume that all the harmonics have the same 
quality factor, Q, and allow the power in each harmonic to be a 
free parameter. We then generate a QPO light curve, Lgpo, using 
these narrow Lorenzians as input to the Timmer & Koenig (1995) 
algorithm, and add this to the light curve already created for the 
broadband noise. 

We show an example of the final predicted PSD in Figure |9] 
using the fiducial model parameters with ro = 50 (black), 20 (red) 
and 10 (green). For clarity, we have set the normalization of all 
the harmonics other than the fundamental to 0, set the width of the 
QPO using the model prediction of ogpo/ fQPO, and set its rms 
power oc l/vo- These PSD show all the main features seen in the 



© 2010 RAS, MNRAS 0Q0.[Tlfl3l 



8 A. Ingram & Chris Done 



-f- 



mi 



-+' 



0.02 



0.04 0.06 
L (arbitrary units) 



0.08 



o 




+ 

+ 


oo 




+ 


s 




+ 


g 




+ 


cs O 




+ 










+ 












■ ++ 











10 



average f (Hz) 



Figure 8. Precession frequency plotted against luminosity, where both are 
calculated at 4s intervals using the fiducial model parameters (gray points). 
After binning (red crosses), we clearly see a linear relationship between the 
two quantities. This relationship has recently been discovered in data from 
the 1998 rise to outburst of XTE J1550-564 (Heil et al in print), demonstrat- 
ing the substantial predictive power of this model. 




f(Hz) 

Figure 9. The full PSD calculated using the fiducial model parameters with 
To = 50 (black), 20 (red) and 10 (green). Here, the QPO is represented by 
a Lorentzian centred at the precession frequency with the width set by the 
rm.s variance in precession frequency (see text). 



data during spectral transitions of BHB (Gierlinski, Nikolajuk & 
Czemy 2008). 



7 FITTING TO XTE J1550-564 

We now have a model which can produce both the broad band 
noise and QPO self consistently from the same changing geometry 
as required by the corresponding evolution of the energy spectra, 
namely a changing outer radius of the hot flow (set by the changing 
inner radius of the thin disc). This model gives the major correla- 
tion between the low frequency break and low frequency QPO. The 
propagating fluctuations naturally give the rms-flux relation, while 
the QPO model also gives a framework in which to understand the 
increase in coherence of the QPO with frequency, and even the 
short timescale correlation of the QPO frequency with flux. Plainly, 
the next step is to use this model to fit real data. 

We use RXTE data from the 1998 outburst of XTE J1550-564 
(Remillard et al 2002; Sobczak et al 2000; Rao et al 2010; Wil- 
son & Done 2001). We look at 5 specific observations with obser- 
vational IDs: 30188-06-03-00, 30188-06-01-00, 30188-06-01-03, 
30188-06-05-00 and 30188-06-11-00; hereafter observations 1-5 
respectively. We only consider energy channels 36-71 (correspond- 
ing to 10-20 keV) in order to avoid disc contamination. 

In all our results so far we have shown power spectra pro- 
duced by averaging the logarithm of the periodogram as Papadakis 
& Lawrence (1993) show that, for a red noise variability process, 
the error estimate converges to Gaussian with only ~ 20 samples, 
rather than the ~ 50 required for averaging the linear power esti- 
mates. Our simulations average over M = 32 samples of the PSD 
but the number of segments from the real data cannot be chosen as 
they are simply set by the length of the available data. These give 
Al = 41, 26, 13, 14 and 14 respectively, implying that the errors 
on the PSD from the data are only approximately Gaussian. A fur- 
ther complication is that the error distribution of the periodogram 
is dependent on the shape of the underlying power spectrum (e.g. 
Mueller, Madejski & Done 2004; Mueller & Madejski 2009) mean- 
ing that even M = 20 may not be enough to give Gaussian errors in 
this particular case. Nonetheless, we choose to use minimization 
for simplicity, but we then do a posteriore checks on the goodness 
of fit using the rejection probability method of Uttley et al (2002) 
and Markewicz et al (2003). This entails first calculating 
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1 
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(AP,,™(/))2 



(10) 



where Psim{f) is the logarithmically smoothed simulated peri- 
odogram, with error APsim{f), to be compared to A4 raw ob- 
served periodograms, Pobs.k{f)- We then re-simulate many (512) 
light curves and calculate many different values of 
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(11) 



The rejection probability, Prej, is given by the percentile of 
Xdis(sijn) values which are lower than the one Xms value. For 
example, if 128 of the XAisisim) values were smaller than the one 
Xdis value, the rejection probability is 25%. This means there is a 
25% chance that the the observed light curve was generated by a 
different underlying variability process than the model. 

One aspect of using real data is that there are uncertainties as- 
sociated with each point on the light curve. Since we use the loga- 
rithmically averaged PSD we need the power to be positive definite 
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so we cannot do white noise subtraction of the data PSD. Instead 
we add white noise to the simulation, and compare this with the 
total (including noise) power spectrum of the real data (see also 
Utfley et al 2002). 

We incorporate our model for the power spectrum into XSPEC, 
using the local model functionality. Our model outputs the loga- 
rithmically averaged total (white noise included) periodogram as a 
function of Fourier frequency rather than the more familiar flux as 
a function of energy. We then fit this to the similarly logarithmi- 
cally averaged total power spectra from the real data. The model 
parameters are the physical parameters governing the behavior of 
the flow. The most fundamental of these is the relation governing 
the viscous frequency at any radius. We have parametrized this as 
a power law function of the Keplarian frequency at each radius, 
fvisc = Br~"^ fkep{r), but the parameters B and m are unknown, 
so are free in the fits, as is F^ar, the intrinsic level of the MRI fluc- 
tuations in mass accretion rate. Then there are the geometric pa- 
rameters To and ri, where also couples to the emissivity, both in 
terms of the index 7 and the boundary condition (either stressed or 
stress free). The combination of Vo, Vi and the surface density pro- 
file fixes the QPO frequency, but we leave its width (parametrized 
as a quality factor, Q) and normalization as free parameters. We 
also allow the normalization of the second, third and subharmonic 
to be free parameters, but fix these to have the same Q as that of the 
fundamental. The final free parameter is the normalization of the 
simulated white noise level, aic, which we then subtract from both 
data and model in order to plot results. 

The truncated disc model predicts that ro is the major param- 
eter which changes across this transition, where remains fixed. 
Hence we fit all 5 power spectra simultaneously in XSPEC, with Vi 
tied across all 5 datasets. Because of the correlation between ri and 
both parts of the emissivity (see section lF!2] l. we choose to compare 
results with some of these parameters fixed at different values. 

7.1 Fit results 

7.1.1 Stressed inner boundary 

In this section we fit our model to the 5 observations considered 
assuming a stressed inner boundary condition. Figure [Tol displays 
the best fit models and data with the irmer radius and emissivity 
index required to be constant across the 5 observations. These ap- 
pear by eye to be good fits, and this is confirmed by their rejec- 
tion probabilities of 51%, 30%, 26%, 12% and 16% for observa- 
tions 1-5, respectively. We find the inner radius to be — 2.45 
which, although low, is greater than the horizon of the black hole 
(rn ~ 1.87 for our chosen a* — 0.5). The truncation radius moves 
from To — 75.7 — 13.6. However, the fit requires that the power 
law index governing the viscous frequency, m, increases across the 
observations, from m = 0.46 — 1.016. We will discuss the signifi- 
cance of this in section[8l 

7.1.2 Stress free inner boundary 

In the previous section, we recovered an inner radius fairly close to 
the horizon (r^ — 2.45) for fits with a stressed boundary condition. 
This assumption implies a causal connection to material closer in 
which, given the proximity of the emission to the horizon, may 
not be physical. Hence the stress free inner boundary condition 
may be more physically appropriate. We repeat the fits assuming 
this, with results shown in Figure[TT] We recover an even lower in- 
ner radius (r^ = 2.0), and similar values for the truncation radius 
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Figure 10. Results of fits for 7 = 4.15 and b{r) = 1 plotted with the white 
noise subtracted. The truncation radius, ro, and rejection probability, Prej, 
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(To ~ 79.9 — 13.8) and power law index governing the viscous 
frequency (m = 0.44 — 1.30). By eye, these again appear to be ex- 
cellent fits; in fact they look nearly identical to Figure[To]and have 
only a marginally higher value. However, these have rejection 
probabilities of 100%, 100%, 87%, 72% and 68% for observations 
1 — 5 respectively, so they are all statistically unacceptable. A re- 
jection probability of 100% means that all of the 512 simulated 
PSDs had X^is{sim) smaller than the one 'observed' xlis' making 
it extremely unlikely that the observed data is part of the simulated 
distribution. 

This discrepancy between and rejection probability points 
to the complex nature of the statistics of the periodogram. We dis- 
cuss this in more detail in the Appendix, but here we simply note 
that these results do not necessarily mean that there is no good fit 
for the stress free version of the model: they just mean that the 
real best fit parameters do not give the lowest . The previous 
fits with a stressed boundary condition also may not be the best 
fit in terms of rejection probability, the difference is that with the 
stressed boundary condition the minimum solution was also a 
statistically acceptable fit on rejection probability. We will explore 
this in more detail in a future paper (Ingram & Done, in prepara- 
tion) 

However, as the fits appear to be good by eye, it seems very 
likely that there should be a set of parameters very close to the ones 
that minimize that give much lower rejection probabilities. 



f(Hz) 

Figure 11. Results of fits for 7 = 5.06 and stress free inner boundary 
plotted with the white noise subtracted. The truncation radius. To, and re- 



8 DISCUSSION 

Although the statistical complexities of the periodogram currently 
limit the model fitting to data, the parameters we derive are still 
most likely very close to the actual best fit parameters, and all four 
fits in the previous section came out with similar results for Vo and 
the viscous frequency radial dependence. While the change in ro is 
predicted from the truncated disc/hot inner flow geometry proposed 
to model the correlated spectral changes, the change in viscous fre- 
quency radial dependence is not. 

There are two (at least) scenarios in which this could occur. 
Firstly, this could indicate that the viscous frequency is well ap- 
proximated by a power law but that the form of this changes as the 
truncation radius moves inwards (as assumed in Ingram & Done 
2010 for the neutron star power spectra). There is currently no the- 
oretical model in which this occurs. A more attractive (to us) pos- 
sibility is that a power law is a poor approximation for the viscous 
frequency and, as the truncation radius moves, the power law that 
best approximates the true shape of the viscous frequency with ra- 
dius changes. This would imply that there is a stationary function 
fvisc{r) (and by extension Vr{r)) which would allow this model 
to fit all 5 observations with only a change in ro. We note that ex- 
actly this behavior is predicted in global simulations of an ADAF 
as the flow has to accelerate strongly to pass through a sonic point 
(e.g. Narayan, Kato & Honma 1997; Gammie & Popham 1998). 
These models give an in-fall velocity f,.(r) which is well approxi- 
mated by a quadratic in log space (see Fig 1 in Gammie & Popham 
1998). This will also allow the model to display a drop off in sur- 
face density in agreement with simulations as discussed in section 
[6] Approximating such a velocity law by a power law of index m 
will result in m increasing as the radial extent of the flow decreases, 
similar to the observed trend. We will explicitly fit such models in 
a later paper (Ingram & Done in preparation). 
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9 CONCLUSIONS 

The truncated disc/hot inner flow model designed to describe the 
spectral evolution of BHB can also give a self consistent geometry 
in which to model the correlated evolution of the power spectrum. 
Propagating fluctuations through a hot flow which extends from 
and outer to inner radius, To — ri, can produce the band limited 
noise characteristic of the continuum power spectrum, as well as 
producing the rms-flux relation (L97; KOI; AU06). Lense Thirring 
precession of this same hot flow can produce the QPO, with fre- 
quency set by the same parameters of Vo and Vi, together with the 
surface density of the flow (IDF09). The surface density is itself 
given self consistently by mass conservation from the propagat- 
ing fluctuations. This predicts that the surface density fluctuates, 
so predicts that the QPO frequency will vary on short timescales 
(i.e. that it is a quasi rather than true period). These fluctuations 
set an upper limit to the coherence of the QPO, and this increases 
(i.e. width decreases) as Vo decreases. This is due to the decrease in 
fluctuation power due to the smaller range of radii from which to 
pick up variability. All these features are well known properties of 
the data (e.g. Remillard & McClintock 2006; DGK07): this model 
gives the first quantitative description of their origin. The fluctua- 
tions also predict that the flux and QPO frequency are correlated 
on short timescales, as a perturbation in the surface density at large 
radii leads to a longer QPO frequency but has little effect on the lu- 
minosity. As this propagates down, it weights the mass distribution 
to smaller radii, increasing the QPO frequency but also increasingly 
contributing to the luminosity due to the centrally peaked emissiv- 
ity. This behavior has also recently been observed (Hell et al 201 1). 

The model also gives a framework in which to interpret some 
otherwise very puzzling aspects of the energy dependence of the 
variability seen in BHB. The extended emission region can be in- 
homogeneous, with different parts of the flow producing a different 
spectrum. The outermost parts of the flow are closest to the cool 
disc, so will intercept more seed photons and have a softer spec- 
trum than that produced in the more photon- starved inner part of the 
flow (Kawabata & Mineshige 2010; Makishima et al 2008; Taka- 
hashi et al 2008). This implies that a larger fraction of the lower 
energy Compton scattered photons come from larger radii in the 
flow than the higher energy ones. The higher frequency variability 
is preferentially produced at the smallest radii, where the spectrum 
is hardest. The flow at these small radii is also furthest from the cool 
disc, so has little reflection spectrum superimposed on the Compton 
continuum. Thus the model predicts that the fastest variability has 
the hardest spectrum and smallest reflected fraction, while slower 
variability has a softer spectrum and larger reflected fraction. This 
trend is also observed in the data (Revnivtsev et al 1999), and is 
very difficult to interpret in any other geometric picture as the in- 
ner disc edge cannot change in radius on even the longest timescale 
(few seconds) over which this relation is seen. 

Similarly, the propagating fluctuations model means that a 
fluctuation starts at larger radii and then accretes down to smaller 
radii. Thus the fluctuation first affects the region producing a softer 
spectrum, then propagates down to smaller radii which produce the 
harder spectrum, so the hard band lags the soft band. The size of 
this lag depends on the frequency of fluctuations considered. Slow 
fluctuations (low frequencies) are produced at the outermost ra- 
dius, so have the longest propagation time down to the innermost 
radius. High frequencies are produced only close to the inner ra- 
dius, so only have a short distance to travel and hence have shorter 
lags. This gives rise to the frequency dependent time lags seen in 
the spectrum (Miyamoto & Kitamoto 1989; Revnivtsev et al 2001; 



KOI; AU06). In future work, we will quantitatively test these ideas 
against real data by incorporating a radial dependence of the emit- 
ted spectrum in our model, as well as putting in the more complex 
form for the radial dependence of the viscous frequency predicted 
from global ADAF models and use improved fit statistics (Ingram 
& Done, in preparation). We will also calculate the light curve 
shape predicted by Lense-Thirring precession of a hot inner flow, 
where the Comptonised emission is modulated by the difference 
in projected area of the (translucent) hot flow, self occultation and 
variation in seed photons from the different projected area of the 
disc (Ingram, Done & Zycki in prep). This predicts non-sinusoidal 
variability in the light curve i.e. this predicts the harmonic structure 
of the QPO. 

While the many successes of the model are clearly evident, it 
is also clear that it is still far from complete. The most obvious out- 
standing issues are of the interaction of the hot flow with the trun- 
cated disc. The mechanism by which the cool disc truncates is not 
well established, though evaporation powered by thermal conduc- 
tion between the two different temperature fluids almost certainly 
plays some role in this (Liu et al 1997; Rozanska & Czerny 2000; 
Mayer & Pringle 2007). Whatever the mechanism, it seems physi- 
cally unlikely that this will give a smooth transition between a cool 
thin disc and the hot flow. Any inhomogeneities will probably also 
be amplified by the difference in velocity between the disc and flow 
(discs are close to Keplarian, while the hot flow is strongly sub Ke- 
plarian) so there will be a shearing turbulent layer formed between 
them. Recent results show that there is variability associated with 
the truncated disc at a few 10s of seconds in the low/hard state of 
the bright BHBs GX339-4 and SWIFT J1753.5-0127 (Wilkinson 
& Uttley 2009), suggesting that there is considerable complexity 
in the disc truncation (see also Chiang et al 2010). Full numerical 
simulations of the MRI in a composite truncated disc/hot inner flow 
geometry are probably required in order to show the effect of these. 
However, such simulations are way beyond current computer ca- 
pabilities. A more tractable issue is the effect of relativity on the 
propagating fluctuations. Near the black hole, light bending and 
time dilation should be important and consequently future versions 
of this model need to take these effects into account. The final goal 
should of course be the creation of a fully relativistic model which 
can produce a Fourier resolved spectrum with both energy and time 
dependence such that we can test it against observations such as the 
PSD, the energy spectrum, the lag spectrum, the cross spectrum etc. 
This is of course very ambitious but it is the only way we can gen- 
uinely achieve a full theoretical understanding of what drives mass 
accretion and emission in BHBs. 
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Figure Al. Raw periodogram points for both data (black) and model (red) 
with case 1 on the top and case 2 on the bottom (see text). As this is the data 
for observation 4, there are 14 black points per frequency bin. Our model 
calculates smoothed periodograms over 32 realisations so there are 32 red 
points per frequency bin. 



APPENDIX A: A CLOSER LOOK AT THE STATISTICS 

Here, we attempt to address some of the statistical problems en- 
countered in section [7] by focusing in on observation 4 as it is a 
particularly striking example. The fit for the stressed version of the 
model (hereafter case 1) is plotted in the 4th panel from the top in 
Figure [TO] and the fit for the stress free version (hereafter case 2) 
is plotted in the 4th panel from top in Figure [TT] By eye, these fits 
both look excellent and are very difficult to tell apart. However, the 
rejection probability is 12% for case 1 and 72% for case 2! This 
is extremely confusing, especially considering that case 2 actually 
has a slightly lower (23.0 compared to 23.8) than case 1! We 
will also discuss the obvious over fitting implied by such low 
values later on in this section. 

In an attempt to understand this discrepancy, we look a lit- 
tle closer at these two examples. It turns out that Xms for case 2 
is larger than in case 1 (59207 for case 2; 53074 for case 1). Al- 
though this sheds a little light on the mystery, it doesn't explain 
such a huge descrepancy in rejection propbability. However, the 
xiis(sim) values are all much smaller for case 2 meaning that 
the simulated values agree with each other much better than they 
do with the observed values, making it extremely unlikely that the 
observed points lie on the distribution of simulated points. In an 
attempt to see this difference, we plot all of the raw periodogram 
points for both cases in Figure lAll The black and red points are ob- 
served and simulated respectively with case 1 on the top and case 2 
on the bottom. We see that the two plots are indistinguishable. 

In Figure lA2l we have looked at 3 specific frequencies: O.lHz 
(black), 9.5Hz (green) and 73Hz (blue) and created a histogram to 
see the distribution of periodogram values. Note, the three distri- 
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we will consider statistical methods more general than but less 
computationally intensive than minimising the rejection probabil- 
ity. Hopefully this will allows us to overcome these statistical prob- 
lems. 
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Figure A2. Histogram showing the distribution of periodogram points at 
/ = O.lHz (black), 9.5Hz (green) and 73Hz (blue) for case 1 (top) and 
case 2 (bottom). The three distributions have been normalised to have the 
same mean for clarity. We can see that these distributions have departures 
from a Gaussian (red line), especially in case 2. 

butions have been normalised to have a mean of and a standard 
deviation of 1 for clarity. Again, case 1 is on the top. Although there 
are departures from Gaussian behaviour in all of these distributions, 
it is clear that the distributions for case 2 are far more centrally 
peaked than those of case 1. This is consistent with ~XAis{sim) be- 
ing smaller in general for case 2. There are also hints of baises 
which may account for the apparent over fitting implied by such 
small values. 

In our analysis, we calculated the logarithmically smoothed 
periodogram over M — 32 realisations in order to calculate and 
minimise x^- We did this because Papadakis & Lawrence (1993) 
show that the logarithmically smoothed periodogram can have a 
Gaussian distribution if it is smoothed over more than AI — 20 
realisations, in comparision to the smoothed periodogram (i.e. no 
logarithms) where you would need more than AI — 50 realisations. 
This implies that our logarithmically smoothed periodograms have 
a Gaussian error distribution, therefore making x^ an appropriate 
fit statistic. However, there are two problems with this, both aris- 
ing because the error distribution depends on the underlying power 
spectrum. Papadakis & Lawrence (1993) were looking at red noise 
which is a significantly different variability process than that con- 
sidered here. Therefore, just because 20 realisations are enough for 
them, doesn't mean it is enough for us. More significantly, the same 
authors show that the logarithmically smoothed periodogram may 
give a baised estimate of a power spectrum with a high curvature. 
Since, we are simulating QPOs, we are in this very regeim which 
could explain the hint of bais displayed in Figure lA2l 

For these reasons, in a future paper (Ingram & Done in prep) 
we will use smoothed periodograms (no logarithms) and attempt 
to find a way around the resulting non-Gaussian errors. To do this. 
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